%拟合形如y=a+bx^2的函数 
%采样点 
x=[19 25 31 38 44]'; 
y=[19 32.3 49 73.3 97.8]'; 
r=[ones(5,1),x.^2]; % 系数矩阵
% 求解系数
ab=lsqlin(r,y);
%设置连续点
x0=19:0.1:44; 
y0=ab(1)+ab(2)*x0.^2; %利用求解出的系数拟合y值
plot(x,y,'o',x0,y0,'r')

y=[3.1 3.27 3.81 4.5 5.18 6 7.05 8.56 9.69 11.25 13.17];
x=0:0.1:1;
x0=[0 0 0 ];
[a,resnorm]=lsqcurvefit(@fit,x0,x,y);
